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ABSTRACT 

The chemical composition of exoplanet host stars is an important factor in understanding the 
formation and characteristics of their orbiting planets. The best example of this to date is the 
planet-metallicity correlation. Other proposed correlations are thus far less robust, in part due to 
uncertainty in the chemical history of stars pre- and post-planet formation. Binary host stars of 
similar type present an opportunity to isolate the effects of planets on host star abundances. Here 
we present a differential elemental abundance analysis of the XO-2 stellar binary, in which both 
G9 stars host giant planets, one of which is transiting. Building on our previous work, we report 
16 elemental abundances and compare the A(XO-2N-XO-S) values to elemental condensation 
temperatures. The A(N-S) values and slopes with condensation temperature resulting from four 
different pairs of stellar parameters are compared to explore the effects of changing the relative 
temperature and gravity of the stars. We find that most of the abundance differences between the 
stars depend on the chosen stellar parameters, but that Fe, Si, and potentially Ni are consistently 
enhanced in XO-2N regardless of the chosen stellar parameters. This study emphasizes the power 
of binary host star abundance analysis for probing the effects of giant planet formation, but also 
illustrates the potentially large uncertainties in abundance differences and slopes induced by 
changes in stellar temperature and gravity. 

Subject headings: planets and satellites: formation — planets and satellites: individual (XO-2) — stars: 
abundances — stars: atmospheres 
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1. Introduction 

The first indication that stellar composition 
plays a role in, and may be affected by, the for¬ 
mation of planets was the enhanced metallicity 

(parameterized as [Fe/Hj]]) of solar-type exoplanet 
host stars (Gonzalez 1997). While initially sug¬ 
gested as a sign of accretion of hydrogen-depleted 
material onto the star, the planet-metallicity cor¬ 
relation is now established as a mostly primordial 
effect (e.g., Santos et al. 2004; Fischer & Valenti 

2005). Recently, the correlation was shown to 
be dependent on planet mass (e.g., Sousa et al. 

2008; Ghezzi et al. 2010), being the weakest 
for terrestrial-sized planets (e.g., Buchhave et al. 

2014). If host star metallicity is considered a proxy 
for heavy element concentration in the protoplan- 


1 [X/H] =Iog(N x /N H ) - log (N x /N H ) so i ar 
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etary disk, this planet mass dependence seems to 
support the core accretion process of planet forma¬ 
tion, as a larger concentration of heavy elements 
is needed to form giant planet cores (~10 Me ; 
Ida & Lin 2004) rapidly enough to allow time for 
gaseous envelope accretion (>100 Me; Pollack et 
al. 1996) prior to gas disk dissipation after a few 
million years (e.g., Wyatt et al. 2008). 

Is it still possible that planet formation may al¬ 
ter host star composition? Melendez et al. (2009; 
M09) suggest the Sun is deficient by ~20% in 
refractory elements, with condensation tempera¬ 
tures T c > 900 K, relative to volatile elements 
when compared to “solar twins” without detected 
planets. M09, as well as studies finding simi¬ 
lar results (e.g., Ramirez et al. 2010; Gonza¬ 
lez et al. 2010), interpret the negative [X/H] 
or [X/Fe] vs T c slopes of refractory elements as 
a star “missing” rock-forming material that has 
been sequestered in terrestrial planets (Chambers 
2010). However, this interpretation has been ques¬ 
tioned. Other studies have found no increase of 
negative slopes in hotter or more metal-rich stars 
(Ramirez et al. 2014); no difference in slopes be¬ 
tween hosts/non-hosts after correcting for Galac¬ 
tic chemical evolution (Schuler et al. 2011b), cor¬ 
relations between the T c slope and stellar age and 
log g (Adibekyan et al. 2014); and even the oppo¬ 
site of expected trends - positive/negative slopes 
in small/giant planet hosts (Gonzalez Hernandez 
et al. 2011). To complicate the picture, stars 
hosting close-in giant planets may instead accrete 
refractory-rich material that would have formed 
small planets (e.g., Raymond et al. 2011; Schuler 
et al 2011b). 

Thus, if/how planet formation influences stellar 
composition remains uncertain, and is challenging 
to disentangle from the local or global composi¬ 
tion of the disk and the star’s position in/motion 
through the Galaxy (e.g., Onehag et al. 2014; 
Schuler et al. 2011b). Binary host stars pro¬ 
vide an opportunity for constraining how planet 
formation affects host stars, as they likely expe¬ 
rience similar environments throughout their evo¬ 
lution (e.g., Kratter 2011). Desidera et al. (2004; 
2006) found that the majority of ~50 wide binaries 
(not known to host planets) had A[Fe/H]<0.03, 
suggesting larger metallicity differences are excep¬ 
tional. Gratton et al. (2001) examined six wide 
binaries across a larger range of elements, find¬ 


ing that four had identical abundances. How¬ 
ever, two stars were significantly different, with 
A[Fe/H]=0.053±0.014 and 0.091U0.006 and simi¬ 
lar differences for other refractory elements; the 
authors suggested these differences were due to 
pollution by dusty protoplanetary disk or rocky 
planetary material. Radial velocity (RV) moni¬ 
toring of both “polluted” systems to search for 
planets has been inconclusive due to changing stel¬ 
lar activity cycles. Within binary systems known 
to host planets, the results are also ambiguous: 
There is not yet consensus on the significance of 
differences between 16 Cyg A and B (the latter 
hosts a ~2.4 Mj planet; Schuler et al. 2011a; 
Ramirez et al. 2011; Tucci Maia et al. 2014), and 
similar studies indicate no abundance differences 
between stars in a single-host (HAT-P-1, Liu et 
al. 2014) and dual-host (HD 20782/1, Mack et al. 
2014) systems. 

Recently, RV monitoring of the southern com¬ 
ponent (XO-2S; 0.98±0.05 M 0 ; Desidera et al. 
2014) of the ~4600 AU binary XO-2 revealed 
two giant (0.26 Mj, 1.37 Mj), slightly eccen¬ 
tric (e ~0.18, 0.15) exoplanets at 0.13 and 0.48 
AU (Desidera et al. 2014). XO-2N (0.98A0.02 
M 0 ; Burke et al. 2007) was already known to 
host a transiting hot Jupiter (0.62 Mj, e=0, a= 
0.04 AU). This system is one of only four known 
dual-planet-hosting binaries (XO-2, HD 20782/1, 
Kepler-132, and WASP-94). Here we expand the 
analysis of Teske et al. (2013) to investigate 
whether there are abundance differences between 
the stars, how such differences depend on derived 
stellar parameters, and how they may relate to the 
formation of different types of planets. 

2. Observations and Data Analysis 

The observations and data reduction for XO- 
2N and -2S are detailed in Teske et al. (2013). 
The R ~60,000, S/N~170-230 data were acquired 
in February 2012 with the High Dispersion Spec¬ 
trograph at the 8.2 m Subaru Telescope, and re¬ 
duced using standard techniques with the IRAF0 
software package. 


2 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni¬ 
versities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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2.1. Original Stellar Parameters 

The stellar parameter analysis and uncertainty 
calculations are also explained in Teske et al. 
(2013) and listed in Table [T| The abundances 
were determined using the local thermodynamic 
equilibrium (LTE) spectral analysis code MOOG 
(Sneden 1973), with model atmospheres interpo¬ 
lated from the Kurucz ATLAS9 NOVER gridtH. 
Equivalent widths were measured with ‘splot’ in 
IRAF, and abundances were normalized to solar, 
derived from a reflected-moonlight spectrum taken 
with the same configuration during the same run, 
on a line-by-line basis. In XO-2N, 49 Fe I and 8 
Fe II lines were measured; in XO-2S, 49 Fe I and 
10 Fe II lines were measured; the line properties 
and equivalent widths are provided in Table [2] 

2.2. Alternative Stellar Parameters 

As shown in recent studies of GAIA benchmark 
stars (Jofre et al. 2014) and transiting exoplanet 
host stars (e.g., Torres et al. 2012), different 
analysis methods lead to different stellar param¬ 
eters and thus different abundances. To investi¬ 
gate this effect thoroughly, particularly the differ¬ 
ent T e ffS, we explored several alternative methods, 
used widely in other exoplanet host star studies, 
for deriving stellar parameters for XO-2N and -2S. 

Alternative Parameter Set 1 First, we re¬ 
peated our analysis relative to solar, using by-hand 
IRAF measurements of the lines in Ghezzi et al. 
(2010). This list includes a smaller number of Fe 
I and II lines (in our case 23 Fe I for N/S, and 
9/10 Fe II for N/S), but all have lab-measured 
and thus more reliable gf-values. Kurucz AT- 
LAS9 ODFNEW grids were used as input mod¬ 
els, though the difference between NOVER and 
ODFNEW are insignificant. The resulting param¬ 
eters are listed as “Alt Param 1” in Table [U 

Alternative Parameter Set 2 Second, we 
adopted the full analysis approach of Ghezzi et al. 
(2010) and derived parameters based on absolute 
abundances, not relative to solar. The Ghezzi et 
al. (2010) line list, with lab-measured gf-values, 
was compiled specifically for an absolute analy¬ 
sis. Again, the measurements were done by-hand 
in IRAF and the stellar models were ODFNEW, 
using the same EWs as “Alt Param 1”. The re- 


3 See http://kurucz.harvard.edu/grids.html 


suiting parameters, are listed as “Alt Param 2” in 
Table Q] 

Alternative Parameter Set 3 Third, we 
used ODFNEW models with a different and longer 
line list of 72/67 Fe I and 9/8 Fe II lines for N/S, 
compiled from Sousa et al. (2008). The equiva¬ 
lent widths were measured with ARES (Sousa et 
al. 2007); lines >0.30 dex from the average in the 
Sun were removed and an interative 2<r clipping 
was applied to the stars during the parameter de¬ 
termination. The resulting parameters are listed 
as “Alt Param 3” in Table [T| 

2.3. Parameter Comparison 

In almost every case, within ler total errors all 
derived parameter pairs for XO-2N and -2S over¬ 
lap, indicating the stars are similar, potentially 
identical, but could also differ by as much as ~200 
K in T e ff, ~0.25 dex in log g , ~0.10 in [Fe/H], 
and ^0.07 km/s in £. In the first three param¬ 
eter sets, T e g of XO-2N < XO-2S, which agrees 
with the original XO-2Nb discovery paper (Burke 
et al. 2007), T e ffS calculated from the IRFM pho¬ 
tometric calibration of Casagrande et al. (2010) 
using colors from Burke et al. (2007), and the 
recent analysis of HARPS-N spectra by Damasso 
et al. (2015). The original and alternative pa¬ 
rameter set 3 both find log g^ > log gs, while 
alternative parameters 1 and 2 find the opposite; 
however, given the log g errors, the magnitude of 
Alog g between the two stars is uncertain (and 
could be zero). Interestingly, a higher [Fe/H] in 
XO-2N is consistent across all derived parameters; 
this also concurs with Damasso et al. (2015). De¬ 
spite the possibility of the stars having the same 
[Fe/H] within errors, the fact that all four parame¬ 
ter sets - found using different lines measurements 
and analyses techniques - find [Fe/H]v >[Fe/H]s 
points towards a real metallicity difference. We 
return to this point in §3. 

2.4. Other Elemental Abundances 

All abundances were derived directly from 
equivalent width (EW) measurements performed 
with the ‘splot’ task in IRAF, except [O/H], which 
was derived through spectral synthesis in MOOG. 
Abundance determinations of Fe, Ni, C, and O 
are detailed in Teske et al. (2013). For elements 
new to this work, line lists (Table [5]) were com- 
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piled from Melendez et al. (2014) and Schuler et 
al. (2011b); atomic parameters are from numer¬ 
ous sources and compiled in VALD (Piskunov et 
al. 1995). All A(N-S) abundance differences are 
listed in Tabic [3] 

The only clement besides Fe for which we mea¬ 
sured two ionization states was titanium. We re¬ 
moved la outliers based on the mean and stan¬ 
dard deviation of the line-by-line solar-normalized 
abundances. This procedure changed the [Ti I/H] 
and [Ti II/H] values by <0.025, but reduced the 
dispersion in the abundances by roughly a factor 
of two. 

Bergemann (2011) found the discrepancy be¬ 
tween solar Ti I and II is solved by including non- 
LTE line formation, and suggested that Ti II may 
be a better [Ti/H] indicator in solar-type stars. 
We tested whether non-LTE corrections to Ti I 
lines in XO-2N and -2S, calculated using statistical 
equilibrium and line formation codes described in 
Bergemann (2011) and MAFAGS-OS models with 
the same stellar parameters derived here, change 
the [Ti I/H] values. We find these corrections do 
not significantly alter [Ti I/H], but as these are 
the first non-LTE corrections calculated for cool, 
metal-rich stars, we provide the A(Ti 1)nlte~ 
A(Ti 1)lte corrections in Table [2] We report a 
total [Ti/H] based on the lines of both species, 
with la outliers removed within each species as 
described above. 

3. Results and Discussion 

The results of our measurements are in Table 
[2j including wavelength, y, log gf , EWs for each 
element for the Sun, XO-2N, and XO-2S. The final 
derived stellar parameters and their ler uncertain¬ 
ties are in Table [U and the (A(X/H)) and a values 
calculated on a line-by-line basis for XO-2S-XO- 
2S are in Table [3] Next we discuss the resulting 
abundances from the four different parameter sets 
derived above. 

3.1. Relative Abundances versus T c 

The advantage of studying binary stars is the 
likelihood that their formation environments were 
similar, so their compositions - assuming like XO- 
2N and -2S they are of similar mass and tempera¬ 
ture - should also be similar. Small differences in 
binary star abundances were previously suggested 


to be related to planet formation/migration (e.g., 
Tucci-Maia et al. 2014). To investigate the XO-2 
system for such abundance differences, we com¬ 
pare one star directly to the other, A(XO-2N - 
XO-2S). In Figure [T] we show the line-by-line (N- 
S) differences and standard deviations of the dif¬ 
ferences (error bars) for our four parameter sets. 

The black asterisks in Figure [0 represent the 
largest possible differences between (N-S) - -204 K 
AT e ff, 0.27 dex Alog g, and 0.11 dex [Fe/H]. Based 
on these differences, XO-2N may have less refrac¬ 
tories than XO-2S, save a few elements - Si, Fe, 
and Ni. These “original” parameter results also 
imply XO-2N is enhanced in volatile elements; the 
slope for T c >500 K elements is negative (—1.66E- 
4 dex/K). 

However, the abundance differences, and thus 
slopes, change with different stellar parameters, as 
suggested in §2. Using a smaller line list, relative- 
to-solar/absolute results in the blue closed/open 
triangles in Figure [l] These parameters corre¬ 
spond to (N-S) differences of -83/-72 K AT e ff, 
-0.09/-0.11 dex Alog g , and 0.09/0.11 dex [Fe/H], 
Interestingly the relative-to-solar and absolute 
analyses produce a very similar pattern, indicat¬ 
ing that XO-2N is enhanced in most refractory 
elements and depleted in volatiles. The red points 
in Figure[l]represent abundances represent param¬ 
eters derived from a longer line list than the “origi¬ 
nal” parameter set, with measurements made with 
ARES, and (N-S) differences of 3 K AT e ff, 0.06 
dex Alog g, and 0.08 dex [Fe/H]. These abun¬ 
dance differences show an even greater enhance¬ 
ment/depiction in XO-2N of refractory/volatile 
elements versus the abundances marked with blue 
triangles. 

The different A(N-S) abundances, and corre¬ 
sponding T c -slopes, resulting from different stel¬ 
lar parameters make it challenging to conclude 
whether any of the enhancements/depletions be¬ 
tween XO-2N and -2S are real, within the bounds 
of measurement errors. To investigate how depen¬ 
dent the abundance differences are on AT e ff and 
Alog g , which change the most between our dif¬ 
ferent analyses, we used the same EWs and three 
test parameter sets that varied only in Alog g. For 
these tests, we adopted AT e ff -130 K (N-S), based 
on the spread of the points in Figure |T] this AT e g 
is between the AT e fjs of the black asterisks (-204 
K) and red circles (3 K). Specifically, the param- 
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eters for XO-2N/S were 5390/5520 K, 0.43/0.33 
dex [Fe/H], and 1.20 km/s for £, with log g vary¬ 
ing between 4.30/4.40, 4.35/4.35, and 4.40/4.30. 

The results of changing Alog g between XO-2N 
and -2S are shown in Figure [2] as blue (-0.1 dex), 
green (0 dex), and orange (0.1 dex) stars.Changing 
Alog g has a minimal effect on the abundance dif¬ 
ferences, except for A(0/H), which is known to 
depend strongly on log g. For most elements, the 
log g tests result in A(N-S) abundances of ~0, 
falling somewhere in between the black, blue, and 
red points from Figure Q] 

However, in all three tests in Figured) the abun¬ 
dances of Si, Fe, and Ni are still above the zero 
difference line. This is remarkably consistent with 
all four parameter sets derived in §2 and plotted 
in Figure [T] which all show Si, Fe, and possibly 
Ni enhancements in XO-2N. Varying many differ¬ 
ent aspects of the analysis and stellar parameters 
does not alter this enhancement pattern, provid¬ 
ing strong evidence for a real abundance difference 
in Si, Fe, and Ni between XO-2N and -2S. 

3.1.1. Possible Causes of Abundance Pattern 

Robinson et al. (2006) and Brugamyer et al. 
(2011) found differences in [Si/Fe] between known 
host and non-host stars; Robinson et al. (2006) 
also found that host stars are enhanced in nickel. 
While these studies compare known hosts and non¬ 
hosts, they point specifically at Si and Ni as im¬ 
portant to giant planet formation. This is inter¬ 
esting given these are two elements we find en¬ 
hanced in XO-2N, the hot Jupiter host. Indeed, 
plotting Brugamyer et al. (2011)’s [Si/H] values 
against planet semi-major axis shows a dearth of 
silicon-poor hosts to close-in (a < 0.1 AU) planets, 
suggesting that silicon may play a role in close-in 
giant planet formation. 

The enrichment of certain elements in XO-2N 
could be due to accretion of material preferentially 
enhanced in these elements. Simulations have 
shown that a significant fraction of giant planets 
migrating in disks as a result of planet-planet scat¬ 
tering cause much, if not all, of rocky disk mate¬ 
rial to be accreted onto the star (Raymond et al. 
2011). The RV evidence for a second object or¬ 
biting XO-2N (Damasso et al. 2015; Knutson et 
al. 2014; Narita et al. 2011) adds weight to the 
planet-planet scattering migration scenario (Rasio 


& Ford 1996; Nagasawa et al. 2008), by which XO- 
2Nb could have migrated inward and subsequently 
had its eccentricity and spin-orbit alignment angle 
damped (Narita et al. 2011). 

XO-2S’s lesser amount of Fe, Si, and Ni may 
instead indicate that these elements are impor¬ 
tant to small planet formation; XO-2S may have 
formed small planets that were ejected. Kaib et al. 
(2013) specifically studied wide binaries, finding 
that ~70% of 1 Mq-1 Mq systems with the sepa¬ 
ration of XO-2N and -2S were unstable and ejected 
a planet within 10 Gyr. Based on dynamical sim¬ 
ulations of wide binary planet systems, Morbidelli 
(2014) found that 73% of planets are ejected, and 
that 60% of wide binaries, versus 41% of isolated 
systems, lost a planet. 

3.2. Comparison to Other Host Star Bi¬ 
nary Studies 

In all parameter cases, a few refractory elements 
(Si, Fe, Ni) stand out as enhanced in the cooler 
XO-2N. This pattern is different from other abun¬ 
dance studies of binary exoplanet hosts. Liu et 
al. (2014) found A[X/H] values consistent with 
zero in the HAT-P-1 G0+F8 binary. Though there 
is not consensus regarding differences between 16 
Cyg A (G1.5) and B (G3), the proposed differences 
(Ramirez et al. 2011; Tucci Maia et al. 2014) ei¬ 
ther show ~0.04 dex uniform depletion with T c , 
or a range of depletions increasing with T c at the 
-0.03 to -0.05 dex level. (We switched A-B values 
to B-A, to stay consistent with cooler-hotter star.) 
Mack et al. (2014), who measured the double-host 
binary HD20781/HD20782 (G9.5/G1.5), show in 
their Figure 3 non-zero differences for O, Al, V, 
Cr, Fe, and Co. Are we seeing an effect in A[X/H] 
trends due to T e g, moving from hotter stars and 
no differences to cooler stars and more significant 
differences? Right now this question is hard to 
answer due to different analysis techniques, line 
lists, and the host star nature of these systems 
(whether or not both stars host planets). In Grat- 
ton et al. (2001)’s Figure 1, plotting A[Fe/H] 
versus AT e ff for stars without confirmed planets, 
there is no clear pattern A [Fe/H] is near zero 
for stars with AT e ff ^10, 90, 275, and 330 K, and 
A[Fe/H] = 0.05 and 0.09 for stars with AT e ff ~120 
and 275 K. Desidera et al. (2006)’s data hint at 
larger A[Fe/H] differences at larger AT e gfS, but 
Desidera et al. (2004)’s data show no real depen- 
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dence of A[Fe/H] on AT e ff. More uniform abun¬ 
dances studies, and a larger sample, are needed to 
understand the temperature dependence of binary 
host star abundance differences and whether they 
are related to planet formation. 

4. Conclusions 

We performed a detailed stellar abundance 
analysis on XO-2N and XO-2S, binary stars of 
similar temperature and mass that host planets of 
different types, to look for effects of planet forma¬ 
tion. Though binary stars are expected to have 
approximately identical abundances, some studies 
report differences between host stars that may be 
due to planet formation. 

When comparing the stars relative to one an¬ 
other, there are hints of abundance differences. 
We test how dependent these differences, and the 
A(X/H) slopes with T c , are on stellar parameters 
using two different line lists, solar-normalized and 
absolute abundance analyses, and measurements 
made by hand in IRAF and automatically with 
ARES. Most of the resulting A(X/H) values vary 
widely, but a few refractory elements - Si, Fe, and 
possibly Ni - consistently appear to be enhanced 
in XO-2N. This pattern is independent of analysis 
method and adopted T e ff and log g value, giving 
confidence to the conclusion that XO-2N and -2S 
have real abundance differences. 

This work adds to the few binary host star anal¬ 
yses, which may prove to be the most useful for 
constraining both the effects of planet formation 
on stellar abundances, and the bulk compositions 
of giant exoplanets. However, it also illustrates 
how sensitive binary host star abundances and T c 
slopes are to derived stellar parameters; such stud¬ 
ies approach the limit of measurement errors and 
should be considered carefully. 

The authors wish to recognize and acknowledge 
the very significant cultural role and reverence 
that the summit of Mauna Kea has always had 
within the indigenous Hawaiian community. We 
are most fortunate to have the opportunity to con¬ 
duct observations from this mountain. We thank 
the referee for their comments and edits that im¬ 
proved the paper. Facilities: Subaru 
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Table 1 

Derived Stellar Parameters and A(N-S) Abundances 


Parameters 

Teff (K) 

error 

log 9 (cgs) 

error 

£ (km s 1 ) 

error 

[Fe/H] 

error 

line list 

relative to solar? 

IRAF/ARES 

Orig Params, XO-2S 

5547 

59 

4.22 

0.24 

1.24 

0.07 

0.28 

0.14 

49 Fe II, 10 Fe II 

yes 

IRAF 

Orig Params, XO-2N 

5343 

78 

4.49 

0.25 

1.22 

0.09 

0.39 

0.14 

49 Fe II, 8 Fe II 

yes 

IRAF 

Alt Params 1, XO-2S 

5523 

49 

4.44 

0.11 

1.35 

0.08 

0.36 

0.04 

23 Fe I, 10 Fe II 

yes 

IRAF 

Alt Params 1, XO-2N 

5440 

69 

4.35 

0.19 

1.29 

0.09 

0.45 

0.06 

23 Fe I, 9 Fe II 

yes 

IRAF 

Alt Params 2, XO-2S 

5447 

47 

4.36 

0.11 

1.08 

0.06 

0.28 

0.04 

23 Fe I, 10 Fe II 

no 

IRAF 

Alt Params 2, XO-2N 

5375 

78 

4.25 

0.20 

1.02 

0.11 

0.39 

0.07 

23 Fe I, 9 Fe II 

no 

IRAF 

Alt Params 3, XO-2S 

5403 

39 

4.40 

0.13 

1.06 

0.07 

0.39 

0.06 

67 Fe I, 8 Fe II 

yes 

ARES 

Alt Params 3, XO-2N 

5406 

32 

4.46 

0.08 

1.11 

0.05 

0.47 

0.05 

72 Fe I, 9 Fe II 

yes 

ARES 









Table 2 

Measured Lines & Equivalent Widths 


Ion 

A 

(A) 

X 

(eV) 

log gf 
(dex) 

EW 0 

(mA) 

XO-2S 
EW (mA) 

XO-2N 
EW (A) 

C I 

5052.17 

7.86 

-1.304 

33.9 

36.4 

33.0 


5380.34 

7.69 

-1.615 

19.4 

22.4 

20.5 

0 I 

6300.30 

0.00 

-9.717 




Na I 

4751.82 

2.10 

-2.078 

13.1 


48.0 


5148.84 

2.10 

-2.044 

12.4 

42.0 

45.8 


6154.23 

2.10 

-1.547 

36.2 

82.9 

89.3 


6160.75 

2.10 

-1.246 

55.4 

103.7 

113.0 

Note 

i.—This 

table is 

available 

in its entirety in 

a machine- 


readable form online. 


Table 3 

Derived A(N-S) Abundances 


A(X/H) 

Orig (N-S) 

Orig cr(N-S) 

AP 1 (N-S) 

AP 1 o-(N-S) 

AP 2 (N-S) 

AP 2 cr(N-S) 

AP 3 (N-S) 

AP 3 cr(N-S) 

A(C/H) 

0.165 

0.007 

-0.013 

0.009 

-0.029 

0.009 

-0.066 

0.013 

A(0/H) 

0.16 


-0.075 


-0.078 


-0.075 


A(Na/H) 

-0.073 

0.006 

0.054 

0.038 

0.069 

0.043 

0.075 

0.019 

A(Mg/H) 

-0.013 

0.025 

0.066 

0.009 

0.079 

0.008 

0.075 

0.020 

A(A1/H) 

-0.042 

0.066 

0.069 

0.010 

0.084 

0.009 

0.087 

0.021 

A(Si/H) 

0.104 

0.032 

0.068 

0.030 

0.075 

0.032 

0.050 

0.029 

A(S/H) 

0.105 

0.061 

-0.057 

0.069 

-0.071 

0.071 

-0.107 

0.079 

A(Ca/H) 

-0.109 

0.050 

0.096 

0.042 

0.114 

0.041 

0.088 

0.041 

A(Ti/H) a 

-0.042 

0.115 

0.022 

0.070 

0.039 

0.084 

0.063 

0.099 

A(V/H) 

-0.086 

0.023 

0.075 

0.027 

0.112 

0.032 

0.163 

0.035 

A(Cr/H) 

-0.113 

0.071 

0.043 

0.043 

0.062 

0.046 

-0.006 

0.051 

A(Mn/H) 

-0.105 

0.033 

0.072 

0.042 

0.090 

0.041 

0.073 

0.028 

A(Fe/H) b 

0.083 

0.067 

0.091 

0.066 

0.109 

0.071 

0.074 

0.077 

A(Ni/H) 

0.061 

0.054 

0.061 

0.056 

0.077 

0.061 

0.075 

0.060 

A(Cu/H) 

0.018 

0.049 

0.038 

0.101 

0.051 

0.112 

0.039 

0.073 

A(Zn/H) 

0.019 

0.002 

-0.003 

0.007 

0.010 

0.008 

-0.051 

0.009 


a Combined Ti I and II average. 
b Combined Fe I and II average. 
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Fig. 1. Top: The A(XO-2N - XO-2S) relative abundances versus T c (Lodders 2003). Black asterisks, blue 
closed/open triangles, and red circles show results from using the “original” parameters, “alternative params 
1”, “alternative params 2”, and “alternative params 3”, respectively. A dotted line shows zero difference. 
All models show a Si, Fe, and possibly Ni enhancement in XO-2N. Bottom: Same as top, showing fits to 
T c >500 K elements; the dashed blue line corresponds to the open triangle points. The red/blue fits have 
positive slopes, whereas the black fit has a negative slope. 
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Fig. 2.- The A(X0-2N - X0-2S) relative abundances versus T c for our test cases, keeping AT e ff for (N-S) 
at -130 K, A[Fe/H] for (N-S) at 0.1 dex, and varying A log g (N-S) from -0.1 (blue stars), to 0 (green stars), 
to +0.1 (orange stars) dex. A dotted line shows zero difference. 
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